import numpy as np
import matplotlib.pyplot as plt

# 定义常数 J
J = 30  # 可以修改 J 的值

# 定义 j 的范围
j_values = np.linspace(0, 10, 500)  # j 从 0 到 10，500 个点

# 计算函数值
function_values = (np.exp(-4*j_values/J))*0.7
# function_values = (J-j_values)/J

# 绘制曲线
plt.figure(figsize=(8, 6))
plt.plot(j_values, function_values, label=r'$e^{-300j/J}$', color='blue')
plt.title(r'Plot of $e^{-300j/J}$ vs $j$')
plt.xlabel(r'$j$')
plt.ylabel(r'$e^{-300j/J}$')
plt.grid(True)
plt.legend()
plt.show()